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Abstract 

A study of hadron pair production mechanism is motivated by the recent observed decays B° — ► 
£)(*)+ K~ K° . One novel phenomenon is threshold enhancement of the kaon pair production. We 
show that these decays in the heavy quark mass limit can be factorized into a generalized form. 
The new non-perturbative quantity is the generalized distribution amplitude which describes how a 
quark-antiquark pair transmits into the hadron pair. A proof of factorization of B° — ► K~ K° 
decays to all-orders is performed by using the soft-collinear effective theory. The phenomenological 
application is discussed in brief. 

PACS numbers: 13.25.Hw, 12.38.Bx 
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I. INTRODUCTION 



Recently, the three body charmful decays B — > D^*'K K^*>° were first observed by the 
BELLE Collaboration 1]. They also found a novel phenomenon that the mass spectra of 
the K~K^° pair peaks near their mass threshold. In [2j, a factorization approach which 
extends the conventional (or say, naive) factorization approach in B-meson two body decays 
[jil was proposed to study B° — > £)W+ K~ K° decays. The transition matrix element of 
B° — > K~ K° is factorized into a B — > D form factor multiplied by a K~K° weak form 
factor. The K~K° form factor is constrained from the experimental data for the time-like 
electromagnetic (EM) kaon form factors. The branching ratios obtained in this approach 
agree with the experiment. The phenomenological success indicates that the generalized 
factorization assumption may be valid as the leading approximation of an asymptotic ex- 
pansion. The purpose of this study is to establish a theoretical foundation of the factorization 
for B-meson multi-body decays, in particular for B — > DKK decays in QCD. 

The factorization approach for the non-leptonic, two body B-meson decays had been 
used for a long time [3]. The basic idea is that a hadronic matrix element of the four-quark 
operator is factorized into a multiplication of two simpler matrix elements which one is rep- 
resented by a form factor and another by a decay constant. An argument based on "color 
transparency" was proposed as a mechanism of the factorization j^]. The QCD factoriza- 
tion method develops the naive factorization into a rigorous and systematic framework jsj]. 
The main ingredient is factorization, the separation of perturbative from non-perturbative 
dynamics, a standard method in perturbative QCD for the hard exclusive processes [6]. The 
naive factorization approach is its lowest order (the strong coupling constant a s ) approx- 
imation. A proof of factorization at two-loop order in B — > Dn decays is given by using 
the methods of momentum regions and Ward identities. The all-orders proof is done in 
the recently developed soft-collinear effective theory |7(. The proof of factorization in the 
soft-collinear effective theory has several advantages: the explicit gauge invariance at the 
classical level can be utilized; the operator language is simpler than the diagrammatic anal- 
ysis; the separation of the perturbative from the non-perturbative dynamics is easily done in 
effective theory because they occur at disparate energy scales. We will use the soft-collinear 
effective theory to give a proof of factorization for B° — > D^*> + K~ K° decays. 

The exclusive, non-leptonic three body B meson decays are usually more complicate than 
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two body case. However, the hadronic physics of B° — > K~ K° decays in the heavy 

quark limit (the heavy quark mass m;, — > oo) may be simpler because the final kaon pair has 

small invariant mass. It is well-known in QCD that the large energy behavior of hadron form 

|~| 

factor satisfies a dimensional counting rule [6|. Consequently, the probability of kaon pair 
production at large invariant mass is suppressed by powers of large energy. This phenomenon 
is called the from factor suppression. Thus the dominant contribution comes from the region 
with small invariant mass. For the production of kaon pair at small invariant mass, it is 
analogous to the production mechanism of two pions in the process 7*7 — > tttt with high 
virtual photons and small invariant mass of pion pair Q]. According to this mechanism, 
K~K° mesons are produced by the hadronization of ud pair emitted from the W-boson 
and can be described by a universal, non-perturbative matrix element which generalizes the 
standard description of a hadron in QCD [6]. The new element is generalized distribution 
amplitude(s) which is the crossed version of of the generalized parton distribution (GPD) 
of hadron |2| . 

In the small invariant mass region, the two kaons are nearly collinear and energetic. 
Thus, the argument of "color transparency" is applicable. Since the ud pair moves fast, time 
dilation effect makes the hadronization of kaon pair cannot occur until the ud pair moves 
far away from the remaining system. The transition ud — > K~K° with a small invariant 
mass is soft and is described by a generalized distribution amplitude of kaon pair. For the 
B — > D transition, the light spectator does not require a hard interaction because B and D 
mesons are both heavy. The energetic and collinear ud pair in a color-singlet configuration 
decouples from the soft gluon interactions. So, the strong interactions between K~K° and 
BD systems occur at short distance and can be systematically calculated in perturbative 
QCD. The above arguments will lead to a factorization form depicted in Fig. 1. Up to 
leading order of ^- (w being the kaon pair invariant mass) 1 , the hadronic matrix element of 
a four-quark operator Qi in the weak effective Hamiltonian is expressed as 

(DKK\Qi\B) = F BD (w 2 ) [ 1 dzH t (z)^ KK (C,z,w 2 ). (1) 

Jo 

Here, F BD {w 2 ) is a B —* D transition form factor at the momentum transfer w 2 , and 
$ (z, (, w 2 ) is a generalized distribution amplitude (GDA) of the kaon pair. H denotes 

1 The mass difference between the heavy quark and heavy meson is neglected. 
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FIG. 1: Factorization of process B — > DKK in the large m B limit. 

a hard scattering kernel which is perturbatively calculable. This factorization formulae is 
a natural generalization of the QCD factorization in B — > Dir decays given in jsj]. Just 
as the B — > Dtt decay, the B — > DKK process provides a clean environment to test the 
factorization in the three body B meson decays. 



II. THE FACTORIZATION IN B — > DKK DECAYS 

For the sake of illustration, our discussion will concentrate on the decay B° — > D + K~K°. 
The extension to B° -> D* + K-K° decay can easily performed. We shall work in the rest 
frame of the B meson. It is convenient to use the light-cone variables p^ = (p + ,p~ ,p±) = 
ri- ■ p-±- + n + ■ p~Y + pj_ with = (2, 0, 0j_) and = (0, 2, 0j_) are two light-like vectors 
which satisfy n 2 + = n 2 _ = 0, and n + ■ n_ = 2. Introduce the total momentum of the kaon 
pair P = pi + p2 where p\^ are momenta of the K~,K° respectively, and invariant mass 
P 2 = w 2 . The momentum P is chosen to be mainly in the "— " direction. Define the 
momentum fraction variable ( = Under the above conventions, we have 

Pb = (m B , m B , 0_l), P d = (m B , rm B , ± ), P = f(r)m B , m B , ± ), 
Pi=f((r]m B , (m B , p ± ), p 2 = f((r}m B , (m B , -p±). (2) 

where r = 7] = and p\ = {((w 2 — m 2 K )/f 2 . We have used the "bar" -notation 

for any longitudinal momentum fraction variable u = 1 — u throughout this paper. One can 
obtain the kinematical constraint on the variable CC — 
We consider the kinematic region 

w 2 < m\. (3) 
4 



This requirement is guaranteed by the form factor suppression at large invariant mass. We 
will give an argument of the suppression later. 

The factorization of B — > decays in the soft-collinear effective theory was studied 
in by using the hybrid position-momentum representation. We will use the position 

space formulation in |l2l . Il3l Il4j ] to demonstrate factorization. We also include discussions 
on the suppression of higher-Fock states and endpoint contributions based on the power 
counting of the soft-collinear effective theory. 

The first step of factorization starts from integrating out the heavy W-boson and the 
hard gluons with virtualities between mw and a renormalization scale /i. The weak effective 
Hamiltonian is obtained as 

Heff = ^K d V cb [C„(aOQo + Cs^Qs] • (4) 

where the four-quark operators are 

Qo = cYil ~ 7s)& - 7sK Qs = cY(l ~ 7s)T A 6 d T „(l - -y s )T A u, (5) 

The Wilson coefficients C , C 8 are given at scale /i ~ m b . 

The low energy effective field theories relevant to our process are heavy quark effective 
theory (HQET) and soft-collinear effective theory. The field degrees of freedom are: the 
heavy quark fields h c vh h h v ; the collinear quark fields £ u , £<j; the collinear gluon field A c ; the 
soft quark and gluon fields q s , A s . For collinear momentum, the off-shellness is p 2 ~ w 2 
rather than Aq CD . The small expansion parameter should be A = ^- (For the c quark mass, 
we assume that it is at the same order of b quark mass m c ~ m&). The Wilson lines are 
indispensable elements in HQET and SCET. The soft Wilson line W s and collinear Wilson 
line W c are very useful to give a gauge-invariant quantities. For example, are invariant 
under the collinear gauge transformations and W}q s are gauge invariant under the soft gauge 
transformations. The definition of Wilson lines and the power counting for the SCET fields 
are given in I 111 . Il2f . 

n 

The coupling of soft gluons to collinear quark is a bit complicated. In |10j, the authors 
introduce auxiliary fields and integrate out the off-shell modes (A, 1, A). The soft Wilson line 
is obtained at the current level. In 14], an approach which is similar to the HQET is used to 



decouple the soft gluon from the collinear quark. 



in leading order of A. We adopt the approach in ^] that the collinear quark field transforms 



'he two approaches give the same results 



as £ — > W s (ii-)£ under the soft gluon interactions. For heavy quark, h v — > W s (v)h v . The 
W s (n) represents that the soft gluons go along the n direction. After these transformations, 
the collinear quark and the heavy quark does not interact with soft gluons. 

The next step is to integrate out the hard mode at order of m& scale. The coupling 
of collinear gluon to the heavy quark leads to off-shellness of order mj which needs to be 
integrated from the effective theory. The tree level matching gives collinear Wilson lines 
which can be represented in a gauge-invariant form as W^. At the loop level matching, the 
factorization formulae will be non-local in position space because the collinear momentum 
component p~ is of the same order as the hard loop momentum. At leading order of A, the 
four-quark operators are matched onto the gauge-invariant operators below 



where k, k' represent or 8 and T k ' = 1 for k'=0, T k ' = T A for k' = 8. Note that k and k' can 
be different because the color-singlet and color-octet currents mix each other by hard gluon 
exchange. The C^srrif,, z, //, /Jf) are position space Wilson coefficients. The //-dependence 
of Cki(srrib, z, /x, fip) cancels the dependence of Ck on the renormalization scale fi. The /ip is 
the factorization scale which separates the hard modes from the matrix element of the kaon 
pair. The appearance of fip starts from two-loop order. In concept, /jp is different from \i. 
Compared to the relevant operators given in [l0|, our formula contains explicit soft Wilson 
line W s (v), W s (v'). It is found that hard-soft momentum region contribution does not cancel 
at two-loop order \^ and can be absorbed in the definition of the B —>■ D transition form 
factor. 

For B° — > D + ir~ decays, the next step is to prove that the soft gluons which attach 
the collinear quarks decouple and then cancel. This cancelation occurs for both the non- 
factorizable and factorizable diagrams for soft gluons. Here, the non-factorizable diagrams 
represents the graphs which attach u, d quarks to b, c quarks. For B° — ► D + K~ K° , the 
case is different. The non-factorizbale soft gluons decouple from the collinear fields and 
vanish due to unitarity of the soft Wilson line W% (n_)W s (nJ) = 1 for color-singlet current 
operator. But, the factorizable soft gluons does not cancel completely because the kaon 
pair system contains the valence ss quarks as well as a lot of see quarks and gluons with 
soft momenta of order of w. The non-cancelation of soft dynamics is similar to B — » D 




x 



[hW c ] (sn + )Wj(n.)T lfl T k 'W s (n.) [w^ u ] (0). 



(6) 
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transitions where the interactions with the spectator quark in B meson are soft dominant. 
So the QCD dynamics of kaon pair is not collinear dominant as the single kaon or pion. 
However, the non- cancelation of soft interactions does not break down the factorization 
because they occur far away from the hard interaction point. The collinear fields in Eq. (JHJ) 
does not receive the non-factorizable soft gluon contributions. So, the factorization formula 
of Eq. © have factorized the hard interaction with virtualities of order of m b and the soft 
interaction with virtualities of order of w. 

Because the color conservation of QCD, a color-octet current can not couple to two color- 
singlet kaons. Although some literatures discuss the possibility of color-octet contribution 
for heavy quarkonium, there is no indication to introduce the color-octet mechanism for 
light hadrons. So the color-octet operator part in Eq. (jUJ) can be set to zero. From the 
above discussions, we can separate the hadronic matrix element (D + K~ K°\Qk\B°) into two 
separate parts as 

(D + K-K°\Q k \B°) = J rf S d(sm b , 2 , / i, / i F )(D + |^,H/ s t (^)r , W s ( V )^| J B°) 

x (K-K°\ [l d W c ] (sn + )T^ [\V^ U ] (0)|0), (7) 

Since £ = ^~f + ifj C} the Dirac spin structure in Eq. (J7|) has only one choice Y = 
in leading order of A. The kaon meson is spin zero, so the axial part of the collinear 
current operator does not contribute owing to parity conservation. The spin matrix ^+7^ 
contribution is possible but it does not appear in leading order of A. Consequently, only the 
quark vector current is left which is different from the case of a single kaon. Similarly, for 
heavy quark Dirac matrix, V — t/l_{1 — 75) in our special case. 

The matrix element for B — > D transition is defined by 

(D+(v%h c v ,WKv'h,W s (v)h b v \B°(v)) = ^^E^{v ■ v')(v + u%, (8) 

The relation between the Isgur-Wise function and the B — > D transition form factor is 
F^ D = ^|==|i £ (y ■ v'). The leading power generalized light-cone distribution amplitude 
for the K~K° pair is defined by the following matrix element jjj: 

(K-( Pl )K°(p 2 )\ [uWc] (sn + )^ [W^ u ] (0)|0) = P, ( l dze izP < s ^^ K - K \z^w 2 ). (9) 

•/ 

where 7 M is dominated by Introducing the Wilson coefficient in momentum space as 

C k (z,ii/m b ,ii F /m b )) = J ds e izP<sn+) C k (sm b , z, fx, fi F ), (10) 



Inserting Eqs. (JB1 H3 HSJ) into Eq. (J7J), we obtain the final factorization form 

(D+K~K°\Q k \B°) =N£(y. v', /i ) C C k (z, fi/m b , fi /fi F , ^ F /m h )^ K ' K \z, (, w 2 , /^).(ll) 

Jo 

where N = y/m^noiv + v') ■ P. Here, we have chosen Ck(z, fi/rrib, [iol /^f, HF/mb) as a 
dimensionless functions. In Eq. (jllj) . we added a scale /io to show explicitly that the Isgur- 
Wise function £(v ■ v', /i ) has a different evolution equation from $ K K °(z, (, w 2 , fi F ). This 
point is not addressed in the previous literatures. 

The Wilson coefficients Ck(z, fi/m b , [i / fi F , fi F /m b ) are infrared finite because they are 
obtained by matching from the full theory onto the low energy effective theory. They do 
not depend on the details of the low energy dynamics. The convolution form of the factor- 
ized form is due to that both the hard coefficient function and the generalized distribution 
amplitude depends on the light-cone momentum fraction z. 

The generalized distribution amplitude $(z, (, w 2 ) provides an important theoretical tool 
to study the production of two hadron pair. It is the time-like version of a generalized parton 
distribution (GPD) of hadron. Since the study of GDA in B decays is only at the start, we 
display some general properties of GDA which is helpful to understand it. 

GDA contains much fruitful physical information. Q(z, (, w 2 ) depends on three variables: 
quark fraction z = ^t, which describes how the current quark shares the total momen- 
tum; hadron fraction £ = -pr, which characterizes the momentum distribution between two 
hadrons; and the invariant mass w 2 . One special feature of GDA is that it is complex in 
general. The imaginary part of $ is due to rescattering effects or resonance contributions. 
The strong phase shift induced by this soft mechanism is neither power nor perturbative (a s ) 
suppressed because the final state interactions between two hadrons occur at low energies. 
Thus it gives an origin of a large strong phase. In this paper, we will not explore this point 
further since only the absolute value of $ is relevant. Another feature of GDA is that it does 
not select all the valence quarks of the kaon pair in the hard scattering at the quark level. 
The additional valence quark pair ss contained in GDA plays no role in the hard scattering. 

The K~K° GDA has only one isospin state, i.e. 1=1. For iso-vector amplitude 
<$> K K °(z,(,w 2 ), the charge conjugation invariance gives 

$(z, C, w 2 ) = $(1 - z, C, w 2 ) = -<S>(z, 1 - C, w 2 ). (12) 

The amplitude is odd under Q < — > 1 — (, so the skewness of the hadron momentum distri- 



S 



bution is described by the £ dependence. The GDA & K K °(z, £, w 2 ) is normalized as: 

fdz <t> K - K °(z,(,w 2 ) = (2C- 1)F™V), (13) 
where ^°(w; 2 ) is the K~K° from factor in the time-like region, which is defined by 

(K-(p 1 )K°(p 2 )\d(0) 1 ,u(0)\0) = (p 1 - P2 ) fl F K - K °(w 2 ), (14) 

The time-like form factor F K K °(w 2 ) needs to be determined from experiment. Eq. (jl3j) 
means that the time-like form factor can be interpreted from a more general concept, namely, 
GDA. 

The GDA Q(z,(,w 2 , p) will depend on the renormalization scale which is at the order 
of 777,6 in our case. Since the scale dependence is only related to the non-local product 
of quark fields, the evolution of $>(z, (,w 2 ,p) is the same as the BLER evolution of the 
pion distribution amplitude [15j. In the limit p — > oo, the GDA <& K K (z,(,w 2 ,p) has the 
asymptotic form [16]] 

$ K - K \z, C, w 2 ) = 6,(1 - z)(2( - l)F K ~ K \ W 2 ). (15) 

Thus the shape of the kaon pair invariant mass spectrum in B° — > D + K~K° in the heavy 
quark limit is completely determined by the time-like weak form factor F K K °(w 2 ). 

In the above derivation of the factorization, we have assumed that the momenta of ud 
quark pair are collinear and the kaon pair are dominated by the small invariant mass region. 
Now, we argue that they are leading power contributions. The endpoint region which one 
quark of ud pair contains most energy while another is soft is suppressed by its phase space 
dz ~ A. The hard coefficients contain only logarithmic dependence of z and does not add 
more power of 1/z. For the transition pion form factor, the consistency of factorization 
requires the suppression of pion distribution amplitude at endpoint. The contribution from 
the large invariant mass of kaon pair is suppressed by 1 /Q 2 with Q the large invariant mass. 
This is the dimensional counting rule for hadron form factor at large Q [6] . Here, we derive 
it from the SCET power counting. 

The power counting is: collinear quark field £ ~ A, kaon meson state \K) ~ A -1 . The 
time-like kaon form factor is described by the matrix element of (0| Q' \K~K°) at large 
Q where the effective operator Q' contains four collinear quarks associated with collinear 
Wilson lines and Dirac matrix elements in leading power. The more fields involved, the 
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higher suppression occurs. We does not investigate the accurate form for the operator Q' . 
The scaling for the (0| Q' \K~K°) is 

(0| Q' \K-K°) ~ A 4 A~ 2 = A 2 . (16) 

Combining with Eq. ()14|l. we derive the well-known result F K K °(Q 2 ) ~ Aq CD /Q 2 . This 
means that the form factor F K K °(w 2 ) at large invariant mass w 2 is proportional to 1/w 2 
and suppressed. This suppression is called the form factor suppression. One important 
phenomenon related to the form factor suppression is that the spectrum of kaon invariant 
mass is enhanced by its threshold region. So, the threshold enhancement of mass spectrum 
is crucial for the consistency of our factorization method. 

The higher Fock-states contribute the subleading power corrections. In the authors 
use explicit calculations to show that 3-particle Fock-state is suppressed by powers of 1/mj, in 
the heavy quark limit. This conclusion can be proved generally from the simple dimensional 
analysis. The m& scalings for the fields and meson state are: quark field £ ~ tt^ 2 , gluon field 

~ ml and meson state \K) ~ rn^ . The scaling dimension of field is determined by its 
ordinary dimension. For kaon meson, we use the normalization (K(p')\K(p)) = 2p°5 3 (p — p') 
to determine its dimension. Adding one gluon field to the effective operators means the 
increase of dimension by 1, the hard coefficients must be suppressed by a power of m& in order 
to match the dimension of m&. Specifically, we consider the gluon field A v = Jq 1 dv vx^G^ 
insertion. The scaling for A v is A v ~ m 2 , so its contribution is 1/m? suppressed compared 
to the leading term. This result is consistent with the calculation in 5(. Note that the above 
dimensional analysis is only applicable for hard interaction. 

Our approach can be considered as the application of the QCD factorization to three- 
body B decays. We call our method as QCD factorization approach. Though it seems that 
our factorization approach is substantially different from the standard pQCD framework 
in their basic ideas are similar. Both are based on the factorization theorem, which 
separates short- distance from long-distance physics in a simple and systematic way. The 
SCET simplifies the proof of factorization. 
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III. PHENOMENOLOGICAL APPLICATION AND DISCUSSIONS 



Now we discuss the phenomenological application of QCD factorization approach into 
B° — > £)(*)+ K~ K° decays. The leading contribution comes from the collinear region where 
the momenta of ud quarks and two kaon are replaced by their largest minus variables, i.e. 
the invariant mass w 2 and transverse momentum are neglected, 

Pl = f(0, Ct|, Ojl), p 2 = f(0, c^§, o±) 

m B _/ 

p n = r(0, z—j=, ± ), p a = r(0, z-j=, 0±). (17) 

The validity of the above collinear approximation needs to be checked by consistency of 
the perturbative result. To the leading power of w/m B , the decay amplitudes for B° — > 
D^+K'K are 

A(B° -> D+K-K ) = ^=V* d V cb a x F^{w 2 ) F*"*V)(2C - l)(m| - m 2 D ), (18) 



and 



A(B° -> D*+K-K°) = ^V: d V cb ai A* D (w 2 ) F K ' K °(w 2 )(2( - l)2m D »e* D « ■ P, (19) 

where Vij are CKM matrix elements, F^ D , A^ D are B — > DW transition from factors defined 
in 0, ai is the Wilson coefficient, and e* D » is the polarization vector of D*. In the above 
equations, we have used the asymptotic form for K~K° generalized distribution amplitude. 

One potential improvement of the QCD factorization approach is that we can calculate 
the matrix elements beyond the tree level. We use the one-loop results for the Wilson 
coefficient a\ given in jf|. The theoretical input parameters are chosen as: a± = 1.05, 
F^ D (0) = Aq D (0) = 0.6. These parameters provide a well fit to the decay modes of 
B° — > D^ + ii~ . The only unknown input is the K~K° weak form factor F K K °(w 2 ). This 
from factor has been constrained via an isospin relation from the experimental data of time- 
like electromagnetic kaon form factors in The obtained F K K °(w 2 ) can be approximated 
as a power-law distribution: \F K K °(w 2 )\ ~ ^4. The discrepancy between this simple model 
and the best fit result in |2| is within 20% level. It should be noted that this determination 
method is not a direct way to extract the K~K° weak vector current form factor. So there 
are still large theoretical uncertainties coming from the form factor. Table |J and Fig. 121 
give the numerical results for branching ratios and kaon pair invariant mass spectra of 
BO D^+K-K decays. 
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1.77 
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FIG. 2: The kaon pair invariant mass spectra of B° — > D + K K° (solid line) and B° — > D* + K K° 
(dashed line). 

Another test of pQCD comes from the ratio of decay rates 

_ T{B^D*+K-K°) ^ ( A (w>) (m|-m^) \ 2 

In the above relation, we have neglected the effects caused by the phase space difference. The 
use of this ratio can reduce the theoretical errors caused by the uncertainties of F K K °(w 2 ). 
In QCD factorization approach, the ratio tk-k° is slightly smaller than 1, while in factoriza- 
tion approach it is about 2 [2j. The difference between the predictions in the two approaches 
lies in the collinear approximation adopted in QCD factorization approach. Whether this 
approximation is reasonable or not is crucial for the validity of applying factorization at the 
realistic energy of m^. 

From Table HI the theoretical predictions of branching ratios are consistent with the ex- 
perimental data. For the ratio tk-k°i it needs further tests. The momentum spectra plotted 
in Fig. |2shows a similar momentum distribution for pseudoscalar and vector D mesons. The 
fraction comes from the range w < 1.5GeV is about 35% which is not sufficient to guaran- 
tee the validity of the collinear approximation. Note that this numerical result is based on 
our insufficient information of the K~K° weak vector current form factor. From physical 
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considerations, the K~K° form factor in the small invariant mass region is likely to be en- 
hanced by the resonance contribution or soft re-scattering effects. The expected momentum 
spectra should be more concentrated in the mass region close to the threshold where the two 
kaons are collinear. This conjecture is reinforced by the experimental measurement that the 
fraction of B~ — > D°K~K° signal events in the invariant mass range w < 1.3GeV is 55%. 
The best fit K~K° form factor in [2| does not satisfy this criterion. If the future experi- 
ment observes that most of the contribution comes from the small invariant mass region, 
such as w < 1.5GeV, it will provide a strong support of our conjecture. Furthermore, we 
suggest to extract the K~K° weak form factor directly from the momentum spectrum of 
B° -> D+K-K . 

One conclusion can be obtained from the QCD factorization is that the power correction 
is proportional to w/rrib- So, the factorization is only applicable for the invariant mass 
smaller than 2GeV. For three body baryonic B decays, the threshold energy of two baryons 
is higher than 2GeV. So, three body baryonic B decays is at the margin of the factorization 
method. Although the K~K° form factor is not known accurately at present, we give a 
crude estimate about the theoretical errors in 

50 _^ D^ + K-R° decays based on power 
counting. The next-to-leading power correction at the amplitude level is proportional to 
which is about 20%. At the decay rate level, the theoretical accuracy within 40% is 
possible to be accessible in QCD factorization method. This accuracy is not as good as that 
in B —>■ Dtt, but it is still important in explaining the experimental data and understanding 
the hadronic physics of three body decays. 

The principle that the hadron pair produced through quark- ant iquark pair can be applied 
to other three body B-meson decays, such as Dim, ttttK etc. For these processes, more 
generalized distribution amplitudes are required. Some studies of using GDAs (or called 
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. The detailed 



two-hadron distribution amplitudes) in B decays have been considered in 
exploration of this subject is also required. 
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